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ABSTRACT 


This study investigates the failure modes, failure strengths, and failure 


criteria of laminated composite plates with stress concentration and subjected 


to bending loading. Graphite/epoxy composites are used for the present study. 


Lamina material properties, such as stiffness and strength, of the composite 


are obtained by experiments. A series of bending tests are conducted for 


laminated, graphite/epoxy composite plates with and without a hole to 


investigate their failure modes and strengths. In addition, finite element 


analyses are performed to compute stress distributions around holes of the 


composite plates subjected to bending loading. Based on the stress 


computation, a couple of failure criteria are examined to predict the failure 


strengths of composite plates with stress concentration. 
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I. INTRODUCTION 


Many of current modern technologies require materials with unusual 
combination of properties that cannot be met by conventional metal alloys, 
ceramics, and polymeric materials. This is especially true for materials that 
are needed for aerospace, underwater, and transportation applications. For 
example, aircraft engineers are increasingly searching for structural materials 
that are light, strong, stiff, and abrasion-, impact-, and corrosion- resistant. 
Rapid increase in the use of advanced composites is primarily due to their 
high strength-to-weight ratio and high stiffness-to-weight ratios, making 
them ideally suited for the high technology construction. Naval structure 
applications of composite materials are typically depending on this weight 
saving when a composite material is used to fulfill stiffness or strength 
requirements. The Navy has used composite materials since at least 1946 
when two 28 foot personnel boats were designed using a laminated plastic. 
Today, composites are used widely in the fleet and include items such as 
submarine hatches and rudders, MK 46 torpedo propellers, minesweeper 


hunters, and patrol crafts [Ref. 1]. 


Failure of engineering materials is undesirable for several reasons; i.e, 
human lives that are put in jeopardy, economic losses, and the interference 
with the availability of products and services. Even though the causes of 
failure and behaviors of materials may be known, prevention of failures is 
difficult to guarantee. The usual causes are improper material selection and 
processing, inadequate design of components, and their misuse. It is the 
responsibility of the engineer to anticipate and plan for possible failure and, 
in the event that failure does occur, to assess its cause and then take 
appropriate preventative measures against future incidents. 

The purpose of this study is to investigate the failure modes, failure 
strengths, and failure criteria of laminated composite plates with stress 
concentration and subjected to bending loading. Both experimental and 
numerical studies are conducted for the study. As an experimental study, 
four-point bending tests are undertaken to obtain lamina material properties 
as well as to examine failure modes and strengths of laminated composite 
plates with or without holes. As a numerical study, finite element analyses 
are performed to determine the stress distributions around the holes. A 
couple of failure criteria are evaluated using both the numerical and 


experimental studies. 


IL LITERATURE SURVEY 


In recent years, there have been significant studies to determine the 
stiffness and strength of laminated composite materials. One of the common 
experiments used to define the responses of composite materials is the 
flexure test, in which the specimen is subjected to either a central load or 
two quarter-point loads. In many cases, results of these tests are computed 
by use of formulas based upon a homogeneous material, although 
bi-directional (0°-90°) specimens are employed. In other words, the geometry 
of the stacking sequence is not taken into account. 

In 1949, Hoff [Ref. 2] presented flexure formulas which, although they 
were based on an elementary strength of material consideration, appeared to 
be quite adequate to describe the bending of unidirectional or bi-directional 
beams. These theoretical equations, in conjunction with flexure experiments, 
can be utilized to calculate intrinsic material properties. Certain difficulties 
are still encountered in the calculations when highly anisotropic composites 


are tested. Pagano [Ref. 3] reviewed Hoff’s equations and expressed them in 


a more convenient form for computational purpose. 








In 1974, Hann and Tsai [Ref. 4] presented graphs to determine the 


laminate stiffness and strength for the [0°/+45°/90°] boron fiber/epoxy and 
declaimed that a set of generalized graphs for all practical laminates of a 
given composite can replace the current limitation of one set of graphs for 
each discrete laminate. 

The first calculation of the effect of a circular hole on the stress 
concentration in a plate was published by Kirsch [Ref. 5] in 1898. He dealt 
with isotropic plates and computed the stress concentration factor of 3. 
Continuous efforts were made by many researchers to work on the stress 
concentration for various cases. To predict the strengths of laminated 
composites containing through the thickness discontinuities, two related 
criteria based on the normal stress distribution were developed. Whitney 
and Nuismer [Ref. 6] considered a circular hole of radius R in an isotropic 
plate of infinite extent with the origin of the x-y axes system at the center of 
the hole. A uniform tensile stress was applied parallel to the y-axis at 
infinity. As the first failure criterion, they assumed that failure occurred 
when the stress over a characteristic distance, d,, away from the 
discontinuity was equal to or greater than the strength of the unnotched 


material. They proposed the second failure criterion stating that failure 








occurred when the average stress over a distance, a,, from the discontinuity 


equaled to the unnotched laminated strength. The criteria included the size 
effect of discontinuity. Later on, Konish and Whitney [Ref. 7] found that 
simply scaling the stress distribution for an isotropic plate containing a hole 
by a ratio of the isotropic and orthotropic stress concentration factor to get an 
approximate stress distribution in an orthotropic plate with a hole might lead 
to significant errors even close to the hole boundary, particularly if the 
integral of the stress distribution was used as a strength criterion in the 
former study. They proposed a simple polynomial which was obtained by 
adding sixth and eighth order terms to the isotropic solution in order to get 
an approximation for the normal stress distribution adjacent to a circular hole 
in an infinite orthotropic plate. 

In 1981, Hoff and Muser calculated the stress concentration factors for 
circular holes in circular composite plates of cylindrical orthotropy when the 
outer edge of the plates were subjected to a uniform uniaxial traction. They 
showed that the stress concentration factor can be plotted in a diagram as a 
function of S, (the compliance in radial direction) and S (two times 
tangential compliance plus shear compliance) if the ratio R of the plate to the 


diameter of the hole was prescribed. They also made a conclusion that the 


stress concentration factors depended primarily on the local arrangement of 
the fiber at the ee of the hole [Ref. 8]. 

Chang and Chang [Ref. 9] presented a progressive damage model to 
assess damage in laminates with arbitrary ply-orientations and to predict 
ultimate tensile strengths of notched laminates. They developed a nonlinear 
finite element program for laminates containing a circular hole. They 
claimed a good agreement was found between the numerical prediction and 


the experimental data. 


Ill. EXPERIMENT 


This section explains the plate bending theory for laminated composites 
to obtain composite material properties from bending tests. The experimental 
apparatus used, the laboratory condition, and the testing procedure are also 
described. 

A. PLATE BENDING THEORY 

Because the flexural test is performed to obtain material properties of a 
graphite/epoxy composite, the plate bending theory for laminated composites 
is reviewed below. 

For the classical plate bending theory, the following assumptions are 
made. It is assumed that the laminate thickness is small compared to its 
lateral dimensions. It is also assumed that there exists a perfect bonding 
between any two laminae. That being so, the laminae are not capable of 
sliding over each other and they have continuous displacements across the 
interface. It is made another important assumption: namely, a line originally 
straight and perpendicular to the laminate midplane remaines so after 


deformation. 


Finally, the so-called Kirchhoff assumption is made which states that 
in-plane displacements are linear functions of the thickness and interlaminar 
shear strains in the x-z and y-z planes are negligible. These assumptions 
reduce the laminate behavior to a two dimensional analysis of the laminate 
midplane. 


The three-dimensional strain-displacement relationships are: 


Based on the Kirchhoff’s assumption, the in-plane displacements are linear 


functions of the thickness coordinate z, 

U = U9 (x,y) + ZF (X,Y), V = VO(XY) + ZF a (X,Y) ceecersesseserseeees (2) 
where u, and v, are in-plane displacements of the midplane. From the 
Kirchhoff’s assumption, Eq. (3) is also obtained. 

exz =F (x,y) + = Os ey2 = Foley) + $= 0 sceunddinveaansaed (3) 

It follows therefore that the vertical displacement of any point does not 


change in the thickness direction. Finally, the three in-plane strain 


components are expressed as 





Now, it is considered a composite made of n stacked layers or plies. Let 


t be the total thickness of the laminated composite plate. Then, the 


constitutive relationship for the k-th layer is 
= [o], = (2), lel, 
where (Q)1 is the transformed reduced stiffness matrix. Substitution of 
Eq.(5) into Eq.(6) yields 
[oly - (Dy,fe°] +10), (4 
Since a laminated composite is in the state of plane stress, there are three 


stress components. The three stress components produce the following stress 


resultants: 





Nx = oxdz 

Ny = [4 oyde 

Nay = 2, xyz 
These stress resultants have the dimensions of force per unit length. These 
resultants give the total force per unit length acting at the midplane. 
Additionally, moments which are produced by the stresses with respect to 


the midplane are also applied at the midplane, and they are written as 
Myx = (5 Oxzdz 
My- {5 oyzdz 
May = f2, txyede 
From above equations, the stress resultants are written as a summation over 


n plies, i.e: 


Ny ky 


kxy 


0 
=| A12 422 A426 |} ey | +] B12 B22 B26 
416 426 466 }} ex, | | 316 326 566 


0 
E | Aj; 432 4j6 |} & By Biz Big | kx 


[N] = [Al[e®] + [51041 
[Mj = [Bf e°)] + (DITA 


— 
Ai Se (Q5j)x@x - 24-1) 











By- 2 P Ope i 1 + seaneee (13) 
and Dy- 31, (Qi), (2) -2%-))I 
Rewriting Eqs.(11) and (12) together gives 
ace ee te cateies (14) 
In the present study, four different laminates, [0°],,, [90°],,, [0°/90°],.. 
and [(+45°/-45)],, subjected to pure bending moments «___ .tudied. Since 
[0°),, and [90°],, laminates have similar stacking properties, they will be put 


into the same case to discuss: 


1. Case I - [0°], and [90°], laminates: 
For both [0°],, and [90°],, laminate composite subjected to pure 


bending moments, Eq.(12) becomes: 
[MI=(DIA] = a asttcttnteeceecenees (15) 
because the matrix [B] will become zero for a symmetric lay-up. 
That is, 
Mx Di, D2 Dig | kx 
| % Fl Bis Bas by | » ee a6) 
Mxy | | Di6 Dr6 Dew || 
Since there is no midplane stretch, 
Ce ee (17) 


Substitution of Eq. (17) into Eq. (16) and inverse of the resultant expression 


gives 


11 





-1 
ex Di; P12 Pig Mx 
ty |--2| Dy2 D22 Do. My 


exy Di6 D265 PE. Mxy 


Because D,, and D,, are equal to zero for unidirectional laminates (0° 


or 90°) and cross-plied laminates (0°/90°), Eq. (18) becomes: 


-] 
ex Di; Diz 0 Mx 
ey beadinee’ 4 D2 D2 0 My 
exy 0 0 D66 Mxy 


A= det Dj D2 0 
0 0 D66 


D*,,= DyDe 
D*,,= -D,2D6 
D*,,= D, Dy 
D*..= D,,D,,-D’,, 
That is, 
ex D22D66 -Y12P66 


eet | DD D11D 
ey Di DaDes 12°66 “11°66 
exy 0 








If subjected to a uniaxial pure bending moment: M,=M,; M,=M,,=0, the 
Strains are expressed as functions of the layers thickness, material 


properties, and the bending moment. In other words, 


ex = “DM ead iveiaseaies (22) 
For the four-point bending test as shown in Fig. 1, 
By a 
Mo 2w 
where : 
p= load (Ib) 


a= length of moment arm (in) 
w= plate width. 


Rewriting Eq.(22) for D,, gives 


z zPa 
D), = =o = 2exw soceccccscceccccaccs (23) 


For [0°],, and [90°], laminates, D,,=2*h**Q,,/3 , D,,=2*h*Q,,/3 , and 





D,,=2*h**Q,,/3 where h is a half of the total plate thickness. As a result: 


= __32pa_ = 3az L. 
alas: aaa 2 (24) 


Since oy = txy = 0. Eq. (6) yields 


oy = 2B 
4h3w 


13 





Figure 1. 4-point bending test 
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2. Case II - [0°/90°],, cross-ply laminate : 
The derivation is similar to that for [0°] and [90°] laminates. 


However, since [0°/90°),, is nv. symmetric, the matrix [B] will not reduce to 
zero and the situation becomes more complicated. 

The stacking sequence and the thickness of each corresponding ply 
are shown in Fig. 2. From Equation (13), the [A], [B], and [D] matrices for 


[0°/90°},, are obtained as when each layer has the same thickness. 


A,,7h*[(Q,,) + (Q,)] ApmAy, 

A,.=2h*(Q,,), A, A,,20 

Ag=2h*(Q,); 

B,,~h’*[(Q,,) +4* (Q,,)/96,  By=h’*[4*(Q,,) + (Q,,)1/96 
B,,=5h’*(Q,,)/96, ) =) = a | rr (26) 
Beo=5h’*(Q,,)/96, 

D,,=h’*[(Q,,) + (Q,)1/3, D,=D,, 

D,.=2h*(Q,,)/3, D,,=D,,=0 

Dg=2h**(Q,,)/3 - 


3. Case III- [+45°/-45*] laminates : 





i 
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Figure 2. The stacking sequences and corresponding 
layer thickness of [0/90] laminate 





For the [(+45°/-45°),], laminated composite plate, the derivation 
procedure is also more complicated than the [0° ] and [90°] laminates. Due 


to symmetry, matrix [B] becomes null and the detail of matrix [D] is 


provided below: 
ex D*y, D*12 Dig | Mx 
ey ae D* 12 D*42 D* 46 My PESErereyrrirr rere (27) 
exy D*16 D*26 D*eg || Mxy 
where : 
D*,,=D,D,-D’2., D* .*-D 2D 66+D Dg, 


D*,,*D,,D.,-D,,D.» D*_=D, ,Deo-D’ 16s 
D*,,2-D,,D.,+D,,D,¢. D**D,,D_,-D’,.,. 
and 
D,,=h’*[(Q,,) + (Qu) +2*(Qy2) +4*(Q,.)1/6, 
D.=D,,, 
D,.=h*[(Qy1) + (Qh) +2*(Qy) -4*(Qee)VG, ssrseecrnseens (28) 
D,.=D,,h’*[(Q,,) - (Q.)1/48, 
Degh’*[(Q,,) + (Qu) -2*(Q,2) 1/6. 
In this thesis, Eq. (24) and the basic equations of material properties are used 
to to compute the elastic modulii E,,, E,, ....etc. The details to calculate 


these elastic modulii are given in the later chapter. 
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B. APPARATUS 

All the bending tests were conducted in a laboratory whose ambient 
temperature is 22.5+2.0°C with an average relative humidity equal to 
41+5%. The MTS 810 Material Test System was used to conduct all the 
bending tests. 

The materials used during the experimental testing were graphite 
fiber/epoxy composites with [0°J,,, [90°],,, [0°/90°],., and [(+45°/-45°),], 
laminates. The experimental samples were cut from different orientation 
laminate plates with nominal dimensions as 12” by 12” square and 0.15” 
thick into small plates with length of 12” and width of 1”. Some samples 
were drilled a quarter inch hole at the center. A fixture was designed to hold 
the composite samples in a simply supported configuration. Each sample 
was positioned on the top of the supporters. A four-point bending test was 
performed as shown in Fig 1. 

Several strain gauges (1/4 inch CEA-13-250UN-350, with gauge factor 
2.120 + 0.5%) were mounted on both sides of the samples. One strain gauge 
was placed at the center of each side of [0°], [(0°/90°),],, and [+45°/-45°),], 
laminae (without hole) for the purpose of obtaining the properties of 


laminated composite. A [90°] laminate sample (without hole) was mounted 


18 





with two strain gauges on each side, one in the longitudinal direction, and 
the other in the transverse direction. 

Additional strain gauges were mounted on each sample with a 1/4” hole 
at the center to measure the strain response of the sample under load more 
accurately. Figures 3 to 6 show the locations of the strain gauges on each 
sample. Tests of these samples were conducted with the MTS material 
testing machine at a crosshead speed of 0.12 in/min ( or 0.305 mm/min). The 
MTS machine provided readings and a force-displacement print-out on the 
x-y plotter (HP 7054A x-y Recorder) for the applied force and the 
measurement of the deflection of the sample. Strain gauge outputs were 
connected to a Measurements Group SB-10 Switch & Balance Unit, and 
readouts, in microstrain, provided by Measurement Group P-3500 Strain 
Indicator. The results were read and recorded manually upto the fracture of 


the sample. 


19 












on the tension side: 


c— gauge # 1 


L/2 L/2 






on the compression side: 


c= gauge #2 


L/2 L/2 











Figure 3. Strain gauge positions for [0], [0/90], and 
[+45/-45] laminae without a hole 


on the compression side: 





Figure 4. Strain gauge locations for [90] laminate 
without a hole 
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on the compression side: 


0.125% 


on the tension side: 


Figure 5. Strain gauge locations for [0] and [0/90] 
laminates with a 1/4” hole at the center 
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on the compression side: 


0.125% 


Figure 6. Strain gauge locations for [+45/-45] 
laminate with a 1/4” hole at the center 








IV. EXPERIMENTAL RESULTS 


This section presents the results as obtained from the experiments. 
While some description of the experimental and calculated results is 
provided here, a more detailed explanation and physical interpretation of the 
results will be given in the later chapter. In order to ensure consistency in 
recording loading responses, samples were placed as uniformly as possible 
with strain gauges in the same positions. The same crosshead speed and the 
same procedures were also employed for each separate test. The results are 


presented graphically and in a tabular form. 


A. SPECIMENS WITFOUT A HOLE 
The four different laminated composite samples without a hole were 
tested to provide a basis to compute the material properties and to compare 


with samples with a hole. Tables 1 to 5 show the loads and strains recorded 


from each experiment for the samples. The load vs. strain relationship is 


plotted in Figures 7 to 14. 





Table 1. results for a [0] laminate without a hole 





*: failure occurred at approximate 1210 Ibs. 


**: see figure 3 for the strain gauge locations. 


Table 2. results for a [90] laminate without a hole 


load{ Ib ) gauge #2 


ee ee ee ee ee 
ee ee ee ee 


*: failure occurred at approximate 28 Ibs. 
**:see figure 4 for the strain gauge locations. 


Table 3. results for [(0/90)],, laminate without a hole 


*: failure occurred at approximate 401 Ibs. 


**: see figure 3 for the strain gauge locations. 





Table 4. results for [+45/-45]-(1) laminate without a hole 


a ees 
aa eee 
107 





*: no breakage occurred. quit recording due to a very large deflection. 


**: see . gure 3 for the strain gauge locations. 


Table 5. results for [+45/-45]-(2) laminate without a hole 





*: no breakage occurred. quit recording due to a very large deflection. 


**: see figure 3 for the strain gauge locations. 
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B. SAMPLES WITH A 1/4-INCH HOLE AT THE CENTER 


Similarly, three different laminated composite samples, [0°], [(0°/90°)],,, 


and [(+45°/-45°)],, with a quarter inch hole at the center of the plate were 
tested. Five strain gauges were mounted with three on the tension side and 
two on the compression side for each sample. The three gauges on the 
tension side were placed with two at a distance of 0.5 in. away from the 
center and another one at 1.25 in. The locations of strain gauges are shown 
in Figures 5 and 6. The results for each sample were recorded in order to 


compare them with the data obtained from the previous section (section A) . 


Table 6. results for a [(+45/-45),,], laminate with a 1/4-inch 
hole at the center 


gauge #1** 
pai im | 0 
| 29.5 | 2997 | 2,683 
| 49s | -sas1_ | 5022 | 5210 | 4.71 
| 695 | -7,834 | -7,686 | 7,008 | 7,066 | 8,180_| 
| 845 | -10,783 | 10,670 | 10,857 | 9,465 | 11,247 _| 
| 90 | 12282 | -12,151_ | 12,320 | 10,619 | 12,748 | 
ee Ee: ee eee Se ee 


*; no breakage occurred. quit recording due to a very large deflection. 


**: see figure 6 for the strain gauge locations. 





Table 7. results for a [0] laminate with a 1/4-inch hole at the center 


gauge # 5 
ee es a re ee 


| 133-996 | -2,070 | 2127 | 2.075 | 1989 
| sz | 2.277 | -2352_ | 2308 | 2353 | 2,254 
| im | 2879 | 2.649 | 2688 | 2,648 | 

| 190.5 | -2882 | 2049 | 2977 | 2,943 | 

| 210.5 | -3,190 | 3,282, | 3,270 | 3,240 | 3,095 
| 230 | 3,509 | 3,575 | 3,542 | 3568 | 3,366 __ | 
| 250 | 3822 | 4133, | 3838 | ates | 3,574 | 
| 290 | 4520 | 4878 | 4480 | 4.954 | 4,057 
| sto | 4853 | 5,227 | 477s | 5320 | #4173 
| 320 | 5.412 | 5.420 | 5551 | 5,500 | #2655 


*; gauge # 5 indicates failure between 310 and 320 Ibs. 





**: see figure 5 for the strain gauge locations. 
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Table 8. results for a [0/90] laminate with a 1/4-inch hole 
at the center 















load (1b) 
| 20s | -ssz_ | 538 | 527s saz | a6 
Patio | -tos7_ | 04a | 072 | 979 
Po] -toz7_ | asoa | sig | ss7_ | a8 
| gos | -2194 | -2133_ | 2026 | 2072 | 1,882 _ | 
| woos | -2,741 | 2670 | 2533 | 2585 | 2,339 
| si 3,567 | 3.462 | 3,286 | 3,342 | 3.014 
| 160s | 4363 | 4.239 | 4006 | 4075 | 3,679 
| iso | -s16s_ | -5038_ | 4.734 | 4.826 | 4,360_ | 


| 220 | -5.998_ | -5.943 | 5.486 | 5,563 | 5,047 
Be a ee 


*; failure occurred at approximate 228 Ibs. 


**: see figure 5 for the strain gauge locations. 


c. CALCULATIONS FOR MATERIAL PROPERTIES 


1. Poisson's ratio v,,: 
By definition, v, is the Poisson's ratio for the transverse strain in the 
j-direction when stressed in the i-direction, which means 
€ 
la reece (29) 
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for a stress applied in the 2-direction. Applying the result of the [90°] 
laminate experiment to Equation(29), the average value of v,, is found to be 


0.02225. 


2. Elastic modulii E,, and E,,, & Poisson's ratio v,,: 
Here, in order to obtain the material properties, E,, and E,, a 


regression analysis using a FORTRAN program was performed to linearize 


the ratios of load vs. strain for [0°] and [90°] laminates. As a result, 
(Plex) ,o = £154.95 x 103 and (Plex) yoo = £7.245 x 103 


Since 


plugging Eq. (24) and Eg. (30) into Eq. (29) results in the following 


relationship: 


After calculations, results for these two elastic modulii are: 








E,," 6.273 x 10° psi ; E,,™ 0.3649 x 10° psi 
From 
Fil 
y —v). 
12 En. 21 
we also get 


V 270.3825. 


3. Shear modulus G.,.: 
A Matlab program was written in order to compute the material 


property, G,,. The experimental result of load vs. strain ratio of the 
[(+45°/-45°),], laminate was compared with the load vs. strain ratio 
calculated from the program using a guessed value of G,,. The Incremental 
Search method was employed to determine the correct guess value for G,, 


iteratively. The shear modulus G,, obtained from this way was 1.045x10* psi. 


D. CALCULATIONS FOR STRESS & STRESS CONCENTRATION 
FACTOR (SCF) 
The failure load for each sample is shown in Table 9. Here, stresses are 


calculated at z= +h or -h (at the top or bottom of the sample) for each 
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sample. From Eq. (25), the failure stress for [0°],,, [90°],,, and [(+45°/-45°),], 


laminae are 


fo 7 1-840x10> psi/per unit length, 


© ofp = 1.26x10> psi/per unit length, 
So0/wo 7 5-08x103 psi/per unit length, 
© 45 /wo 7 61.86x10> psi/per unit length, 


045; = 29.67x10> psi/per unit length. 


where subscripts wo and w indicate “without a hole” and “with a hole”, 
respectively. 


Table 9. dimensions and failure loads of samples 


*: failure due to a very large deflection. 





For [0°/90°],,, the matrices [A], [B], and [D] were calculated from Eq.(26). 
Plugging these data into Eq.(14) and Eq.(7), the failure stress for the 


[0°/90°],, laminate are 
5.0/90-toplwo ~ 213.86x103 psi/ per unit length, 


©0/90-bottomiwo * 12.53x10 psi/ per unit length, 
50/90-top/w ~ 162.13x10° psi/per unit length, 


50/90-bottomiw 7 9-49%10° psi/per unit length. 


Based on these data, the experimental stress concentration factors (SCF) for 
laminated composites are 1.46 for the unidirectional plate and 1.32 for the 


cross-ply plate if the failure occurs due to the largest stress at an edge of the 


hole. 
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Figure 7. Stress- strain diagram for [0] laminate without a hole 
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Figure 8. Stress- strain diagram for {0} laminate with 


a 1/4" hole at the center 
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Figure 9. Stress- strain diagram for [90] laminate without a hole 
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Figure 10. Stress- strain diagram for [0/90] laminate without a hole 
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Figure 11. Stress- strain diagram for [0/90] laminate 
with a 1/4" hole at the center 
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Figure 12. Stress- strain diagram for (+45/-45] laminate without a hole 
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Figure 13. Stress- strain diagram for [+45/-45] laminate without a hole 
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Figure 14. Stress- strain diagram for (+45/-45] laminate 
with a 1/4” hole at the center 
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V. FINITE ELEMENT ANALYSIS 


A. OVERVIEW 

Since the earliest development of the finite element method, a large of 
plate bending elements have been developed and reported in the literature. 
For example, Batoz et al. [ Ref. 10], with the purpose of identifying the most 
effective element for the thin plate analysis, presented an assessment of flat 
triangular plate bending elements with displacement degrees-of-freedom at 
the three corner nodes only. They claimed that the most efficient and reliable 
three-node plate bending elements were based on the discrete Kirchhoff 
theory and a hybrid stress model. Later on, he [Ref. 11] also presented the 
explicit expression of the stiffness matrix of the nine degree-of-freedom 
plate bending triangular element (called discrete Kirchhoff triangle) which 
allowed, as he said, a significant reduction of algebraic operations in the 
evaluation of the stiffness matrices and bending moments. 

In 1983, Fricker [Ref. 12], based on Razzaque’s report, described a new 
three-noded triangular element for plate bending with an extra internal 


‘bubble’ function proposed by Irons for making the Razzaque's element more 
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flexible. He also claimed that his element and that developed by Hansen, 
Bergan, and Syvertsen were the two most accurate triangular thin plate 
elements currently available. Furthermore, Lee and Zhang [Ref.13] 
developed a six-node plate bending element based on a modified 
Hellinger-Reissener principle and the Reissner-Mindlin plate bending theory 
to obtain more accurate and reliable solution. 

In this thesis, a triangular plate bending element, which was developed 
by Tocher{Ref. 15], is used to analyze the laminated composite plates 
subjected to a pure bending load and to determine the deflections, strains, 


and stresses of the plates. 


B. TRIANGULAR PLATE BENDING ELEMENT FORMULATION 

The plate, in general, may have any irregular geometry and loading. 
The derivations are based on assumptions of the small deformation theory, 
described in the earlier chapter. A triangular plate element is employed for 
the following developement. Each nodal displacement of the element 


possesses three components: a displacement in the z direction, w; a rotation 


about the x axis, w, (derivative of w with respect to y); and a rotation about 





the y axis, w, (derivative of w with respect to x), as shown in Figure 15. The 


displacement function, w, is assumed to be 


Wx,y) = 41 +a5x + any tagx? + agxy + agy? +a7x3 + ag(x2y + xy) +4 


err (33) 
where 
(n} = (1 xy x xyy’ x? (x’ytxy’) yy Jo cee (34) 
and 
a 
a2 
[2] ae er eee (35) 
ay 


Here, the constants a,, a, ...... » a, must be determined from the nodal 
conditions 
W(X,y)™qy> W,(X,y)™Go, -W,(X,y)=q, at (x,,y,)=(0,0) 
W(X,Y) "Gyr W,(Y)* ss -W,(XY)™Gg at (Xp,Yq)=(O,Yq) — sascscecessececees (36) 
W(XY)*Gy> W,(X,Y)™ Gg, -W,(X,Y)"qy at (Xs¥5)™(X,,Y5) 
The local y axis is taken the same as the line connecting the node 1 and node 
2 with the origin placed at node 1. The local x axis is taken towards node 3 
as shown in Figure 16. By using Eq.(33), Eq.(36) can be shown in the matrix 


form as 


45 


q2=Q3i-1 


Figure 15. Nodal degrees of freedom of a triangular plate 
in bending 


Figure 16. Local nodal coordinate system 








GOS | ce [x tea (37) 
q9 
where [n’] is a 9x9 matrix (see Appendix). By using Eqs.(33) and (37), 


Eq.(4) can be expressed as 
= [aa SY rr (38) 
where 
_ 0002006 2% oO 
[#]-- 0000020 rw Sy] ersssscssssseee (39) 
000020 0 4&«+y) 0 


and 


H] = [ H|['y* suet (40) 


Finally, the element stiffness matrix in the local (x,y) coordinate system can 


be written as 
[KO] = 4, 21 (6) 21 (41) 


where V indicates the volume of the element and the matrix [Q] is the 
material property matrix for the generally orthotropic lamina. By substituting 


for [H] from Eq.(40), Eq.(41) becomes 


[x7 -[ (IT) al J) go [i] 70) Bala Audbisee (42) 


where t is the thickness of the plate. For a laminated composite plate 
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cd 3_ 43 
frig? 1Oide~ 5 (Ole #4 -#4-1} ~ I 
and the components within the bracket of Eq.(39) won't vary with the 
thickness. The element stiffness matrix in the local (x,y) coordinate can be 


rewritten as 


[x10] =([o"y 1) [great oucras [ory 


000200& 2% oO 
(C]=-|0000020 2% by 
000020 0 &x+y) 0 
The area integrals appearing inside the bracket of Eq.(44) can be calculated 
in terms of the global XY coordinate system as well as in terms of local xy 


system chosen in Figure 16 from the following relations: 


area tdy = 3.3.99 

Si, 9- 57392 

fy, 2 4 x392.02 +y3) 
2. dedy = 57392 

xy. dedy = > se x5 7202 +273) 


1 2 2 
Gy". dedy = 753.9202 +23 +3) 


Appendix shows the detailed saucuiations of the local element stiffness 





matrix. Finally the element stiffness matrix from Eq.(44) in the global 


coordinate system whose XY plane is assumed to be same as the local xy 


plane) can be obtained from 
(KY) = DT Yop aaeeeeeee (52) 
where the transformation matrix [p] is given by 
10 000 000 90 
0 lox mox 0 0 0 00 O 
0 loy moy 0 0 0 00 0 
00 010 000 0 
(PJ=|} 0 0 O O lox max 0 0 OD | aassccereccccccsccceces (53) 
x9 100 0 Oloymoy0 0 0 
00 000 0310 0 
00 000 0O O lox mox 
00 000 0 0 loy moy 


where (15x, Mox) and (I5y, Moy) represent the direction cosines of the lines ox 
and oy respectively. 
After each element stiffness matrix is calculated, it is assembled into 


the global stiffness matrix. The resultant matrix equation becomes 


(Alle) =(FE CC adds dase ssnssaavsnne (54) 
where [F] is the force matrix. 


In order to compute stresses and strains, post-processing is 
accomplished. After obtaining displacements, local strains are extracted 


from the global displacements. To this end, the local coordinate, matrix [n’], 
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direction consines, and the transformation function for each element are 
calculated. Then, the local strain vector is 


le ce ei (55) 


where matrix [H’] is the same as matrix [H] except that x is one-third of 
(x,+x,+x,) and y is one-third of (y,+y,+y,) and z=+h and -h in order to 
compute strains at the top and bottom sides of the plate. In addition, [d] is 
the nodal displacement vector of an element. The local strain is transformed 


into the global strain by using the following equation: 


OS Fae Cs (56) 
where 
2, mex lox * mox 
[7] = mex Py lox * mox 


-2lox * mox 2lox * mox 2. - mex 


At the end, stresses are computed from the strains using Eq. (6). These are 


the global stresses. 


50 


VI. NUMERICAL RESULTS AND DISCUSSION 


A. NUMERICAL RESULTS 

In order to calculate the strain and stress for a laminated composite 
plate with a quarter inch hole at the center and subjected to a pure bending 
moment, a program named MFEM.M is written in Matlab. With an input for 
material properties, plate dimensions, the applied bending moment, and the 
failure load of the material, this program calculates the deflection, strain, and 


failure stress for each element. Another program named DMTRX.M is also 
written in Matlab which is capable to calculate the stiffness matrix (Q} for 
four different kinds of laminae [0°], [90°], [0°/90°],., and [(+45°/-45°),].. 


In this thesis, a laminated plate with a quarter inch hole at the center is 
modeled with 214 nodes and 371 elements, which generates 642 by 642 
stiffness matrix (642 dof) as shown in Figure 17. After calculations, 
normalized stresses of stress concentration and the normalized failure 
stresses v.s the distance away from the circular hole edge for [0°)],, and 
[0°/90°},, are plotted in Figure 18 through 23. As a result, the stress 
concentration factor calculated from the F.E. analyses for the [0°],, laminate 


is 1.47 and that for the [0°/90°],. is 1.27. Normalized failure stresses for the 
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[0°},, laminate, the top [0°] layer and bottom [90°] layer of the [0°/90°],, 


laminate are all close to one at the edge of the hole. 


B. DISCUSSION 
Whitney and Nuismer (Ref. 7] said that when the stress at a 
characteristic distance, d,, away oni the discontinuity would be equal to or 
greater than the strength of the unnotched material, failure would occur. In 
this current work, the constant value, d,, for the [0°],, and [0°/90°],, laminae 
is zero from the results. Thus, the failure criterion of Whitncy and Nuismer 
assumption is not suitable in this kind of graphite fiber/epoxy. Since it is of 
great importance for design purpose to be able to predict the strength of a 
composite under loading conditions prevailing in service, some failure 
criteria like the maximum stress criterion and the Tsai-Hill failure criterion 
are used in order to predict the failure and check the failure criteria. 
1. Maximum Stress Criterion 
This criterion tells that failure will occur when any one of stress 
components is equal to or greater than its corresponding allowable strength. 


Since the stress in the x-direction is the dominant stress in this current work, 


failure will occur when 





o} 2x7 


where X", is the ultimate uniaxial tensile strength in the fiber direction ( in 
this case, x-direction). 
2. Tsai-Hill Criterion 
According to the Tsai-Hill criterion, failure of an orthotropic lamina 


will occur under a general stress state when 


2 2 
pad Td a ad hd Ve 
x2 y2 y2~ 2 


where X, Y, and S are the longitudinal tensile failure strength, the transverse 
tensile failure strength, and the in-plane shear failure strength, respectively. 

Figures 7 through 14 show that for this particular laminated composite 
failure occurred almost simultaneously on both compressive and tensile 
sides. From figures of normalized failure stresses, it is obvious that failure 
occurred at the hole edge when the stresses at the hole edge were equal to 
their corresponding allowable strengths. 

Stress concentration factors obtained from the experiments and 


numerical analyses are nearly the same for both [0°],, and [0°/90%),, 





laminated plates. This result confirm that failure occurs when the maximum 


stress at the edge of a hole is the same as the material strength. 
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Figure 17. x-y plane view of the FEM mesh of the composite 
plate with a circular hole. 
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Figure 18. Stress concentration factor for [0] laminate 
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Figure 19. Stress concentration factor for [0/90] laminate at top layer-[0] 
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Figure 20. Stress concentration factor for (0/90] laminate at 
bottom layer-(90] 
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Figure 21. Stress distributions for [0] laminate 
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Figure 22. Stress distributions for [0/90] laminate at top layer-[0} 
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Figure 23. Stress distributions for [0/90] laminate at bottom layer-[90] 
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VIL CONCLUSIONS 


In this thesis, a four-point bending test is used to investigate the failure 
strain and failure stress for a graphite-fiber/epoxy composite of four different 
laminated plates. From experimental results, composite elastic modulii are 
obtained and the stress concentration factors due to a quarter inch hole 
located at the center of the plate are calculated. Furthermore, a numerical 
analysis using the finite element method is conducted to simulate the plate 
under the same loading condition. By using both experimental and numerical 
results, the characteristic distance d,, which determines failure based on the 
Whitney and Nuismer theory, is examined. Finally, the maximum stress 
criterion and the Tsai-Hill failure criterion is applied to check the failure of 
composite. These failure criteria confirm that failure occurs when the 
dominant stress component of the laminated composite is equal to its 
maximum allowable strength. From this work, it is also concluded that: 

¢ the four-point bending test is a good way to investigate material 

properties although it is a little more complicated to conduct than the 


tensile test. 
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¢ The characteristic distance for the Whitney and Nuismer theory is 
zero for the present composite material. Both the maximum stress 
criterion and Tsai-Hill failure criterion are proper to predict the 


failure for this material. 
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APPENDIX 


The [n’] matrix in Equation (37) is 


100 0 0 0 0 0 0 
0010 0 0 0 0 0 
0-10 0 0 0 0 0 0 
10y 0 0 ys 0 0 y3 

2 

ae 001 0 0 yw, 0 : 3y3 
0-10 0 “y2 0 “Y5 0 


1 x3 y3 x3 X3Y3 yz x3 (x33 +x3y3) y3 
00 1 0 x3 2y3 0 (2xgy3+x3) 3y2 
0-1 0 -2x3 -y3 0 -3x2 -(y3+2x9y3) 0 


The detail calculations of [C]"[D}[C] in Equation (44) are shown as 


00 O 

00 O 

00 O 

20 0 Dj; Dj2 Di¢ 
Ic[Dl=| 0 0 2 Dy2 D2 D6 

02 Oo Dig Dog D 

an ae 16 926 P66 |3,3 

2y 2x 4(x+y) 

0 6y 0 9x3 








0 0 0 
0 0 0 

0 0 0 

2D}; 2Dj2 2D16 

[c)"[p1}=| 2P16 2026 Dee 
2D}2 2D22 2D. 

4; 42 &43 
6yD 19 6yD29 6yD76 


where 
Ay = 2yDj1 + 2xD 12 +40 + 16 
Ag = 2yDj2 + 2xD72 + 4% + yD 76 
A3 = 2yD16 * 2xD 76 + 4% + YD66 
and 


000 O 0 0 0 0 0 
000 0 0 0 0 0 0 
000 0 0 0 0 0 0 


000 4D)] 4D16 4D}2 12xD} Ay 12yD 19 
(Cl? [DIC] a 000 4Di6 4D66 4D6 12xD 36 As 12yD 76 
000 4D)2 4D 6 4D>9 12xD}> A6 12yD99 
0.00 12xD,1 12xD,6 12xD)7 36x2D,1 7 36xyD}> 
000 Ay, As Ag Az Ag Ag 
0.0.0 12yDj> 12yDy6 12yD27 36xyD17 Ag 36y722 


Ag = (4D12 + 8D 6) + (4D) ) + 8D 6 
As = (4D 6 + 8D66)x + (4D 16 + 8D E66 
Ag = (4D 29 + 8D26)x + (4D 1 + 8D26)y 
Aq = (12D1 + 24D 6)x* + (12D) + 24D 6 )xy 


Ag = (4D99 + 16D 6 + 16D6g)x2 + (8D 9 + 16D 16 )xy 
+(16D 6 + 32D 66 xy + (4D 1 + 16D y6 + 16D 66)xy 


Ag = (12D 29 + 24D96)xy + (12D 19 + 24D 6)? 


Thus, the area integrals inside the bracket of Eq.(44) can be obtained 


0000 00 00 0 
0000 00 0 0 0 
0000 00 00 0 
000 4; Ay Ag Aq Aq] A16 
gicipitcias =| 900 Ay A3 As Ag Ay2 447 
A 000 Ag As Ag Ag Aq3 Ag 
000 Az Ag Ag 419 414 A19 
000 Aq} 442 413 414 415 A420 
000 Aj6 417 Aig 419 420 421 
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where : A,=2D,,%:Y2 > A,=2D,,X3y2 >» 


A,=2D,,%,y2 » A,=2D,,%,y; » 
A;=2D,,%,y, » A,=2D,,%,, » 
A,=2D,,x*,y2» A,™2D,,*°3y2» 
A,=2D 2X37 » Ao23D1%'sY2 » 


A,,=(4D,,+8D,,)(,¥,/6+(4D,,+8D,,)[X,y2(¥2+Y3)/6], 
A,,=(4D,,+8D,,)\(x",y,)/6+(4D,,+8D,,)[x,y,(y.+¥3)/6, 
A,3=(4D.+8D3.)(x',Y)/6+(4D,,+8D2)[X,Y2(¥2+Y3)/6], 
A, ,=(12D,,+24D,,)(x°,y,)/12+(12D,,+24D,,)[x2,y,(y,+2y;/24], 
A,;=(4Dz,+16D,,+16D,.)(x?,y,)/12 
+(8D,,+16D,,+16D,,+32D,,)[x’,y,(y.+2y,)/24] 
+(4D, ,+16D,,+16D,,)[x,y,(y,+y,Y3t 73/12], 
Aye 2D 2%sY2(¥2t¥3),  Aqy® 2D 2gXaY2(¥2* Ys)» 
Ajis72D.,%,Y(¥.+¥3))  Ajo=3D,.x",Y2{y.+2Y5), 
Ayo=(12D2.+24D,.)[X'3¥2(¥2*+2y3)/24] 
+(12D,,+24D,,)[x,y,(¥'2t 27st 312], 
Ay:=3D 2% Yo(¥'2* 2st Y's)- 
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